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I.  INTRODUCTION 


For  some  applications,  a  decrease  in  the  amount  of  coolant  required 
for  a  transpiration-cooled  reentry  vehicle  nose  tip  is  desirable  and  potentially 
very  beneficial.  The  work  presented  here  seeks  to  optimize  the  shape  of  a 
reentry  vehicle  nose  tip  in  order  to  minimize  the  total  heat  transferred  (and 
thus  the  transpiration  coolant  requirements)  to  the  nose  tip  throughout  a 
realistic  reentry  trajectory.  The  nose  tip  shape  is  assumed  to  remain 
constant  throughout  the  trajectory. 

The  majority  of  work  devoted  to  the  determination  of  optimum  body 

shapes  in  a  hypersonic  flow  environment  has  been  concerned  with  the  problem 

of  finding  minimum  drag  bodies*.  Recently,  however,  work  has  appeared 

which  considers  the  determination  of  minimum  heat  transfer  body  shapes  in  a 

constant  free  stream  environment.  Laminar  flow  solutions  were  found  by 
2  3 

Belyanin  and  Aihara  .  Belyanin  discussed  the  relative  heat  flux  to  his 

minimum  heat  transfer  solution  bodies,  but  did  not  give  any  details  of  the 

4 

body  shapes.  Aihara's  solution  has  been  shown  by  Hull  not  to  be  a  minimum 

5 

solution.  Perminov  and  Solodkin  obtained  nons lender  body  solutions  for 
both  minimum  total  drag  and  minimum  heat  transfer  for  both  laminar  and 

*Miele,  A.  (ed.),  "Theory  of  Optimum  Aerodynamic  Shapes,  "  Applied  Mathe- 
matics  and  Mechanics,  9,  Academic  Press,  New  York  (1965). 

2 

Belyanin,  N.  M.,  "Determining  the  Form  of  a  Body  Resulting  in  Minimum 
Heat  Flux,  with  Laminar  Flow  in  the  Boundary  Layer,  "  Izv,  AN  SSSR, 
Mekhanika  Zhidkosti  i  Gaza,  2j6),  pp.  37-45  (1967).  English  Translation, 
Fluid  Dynamics,  Consultant's  Bureau,  New  York  (1969). 

3 

Aihara,  Y.,  "Optimum  Body  Geometries  of  Minimum  Heat  Transfer  at 
Hypersonic  Speeds,  "  AIAA  Journal,  Technical  Notes,  £(11),  pp.  2187-2188 
(November  1968). 

4 

Hull,  D.  G.,  "On  Hypersonic  Shapes  of  Minimum  Heat  Transfer,  11  The  Journal 
of  the  Astronautical  Sciences,  XVII(l),  pp,  60-62  (July-August  1969). 

5 

Perminov,  V.  D.,  and  E.E.  Solodkin,  "Axisymmetrical  Bodies  with  Minimal 
Resistance  and  with  Minimal  Flow  Toward  the  Surface  of  the  Body,  with  Dif¬ 
ferent  Characters  of  the  Flow  in  the  Boundary  Layer,  "  Izvestiya  Akademii 
Nauk  SSSR,  Mekhanika  Zhidkosti  i  Gaza,  (2),  pp.  94-102  (March-April  1971). 
English  Translation,  Fluid  Dynamics,  Consultant's  Bureau,  New  York  (1973). 


turbulent  flow.  They  found  the  laminar  and  turbulent  minimum  heat  transfer 
shapes  to  be  very  similar  to  one  another  for  a  given  fineness  ratio.  The  work 
presented  here  extends  this  previous  work  to  the  case  of  a  reentry  trajectory 
environment  and  seeks  to  minimize  the  total  heat  transfer  to  the  nose  tip  for 
the  entire  trajectory. 

Section  II  first  presents  a  brief  derivation  of  the  total  heat  transfer 
relationships  used  in  this  work  and  the  approximate  flow  field  evaluation  pro¬ 
cedures  utilized.  Then,  in  Section  III,  the  mathematical  formulation  of  the 
minimum  heat  transfer  problem  using  a  variational  calculus  approach  is  dis¬ 
cussed.  Calculated  results  are  then  given  which  discuss  both  minimum  and 
maximum  heat  transfer  solutions  and  the  relative  heat  transfer  rates  to  these 
shapes  compared  with  several  other  body  shapes. 


II.  HEAT  TRANSFER  PREDICTION  METHODS 


A.  INTEGRATED  HEAT  TRANSFER  EXPRESSIONS 

When  expressions  for  the  laminar  and  turbulent  convective  heat  transfer 

6  7 

coefficient  given  by  Vaglio- Laurin  '  1  are  used,  and  the  procedure  of  Allen 

g 

and  Eggers  is  followed,  the  rate  of  change  of  integrated  heat  transfer  with 
altitude  dQ/dy  is  given  by 


s 


where  n  =  1  for  laminar  flow  and  n  a  1/4  for  turbulent  flow.  When  the  inte 
g ration  in  nondimensional  surface  distance  s  =  s/r^  is  performed,  this 
expression  becomes 


dQ 

dy 


-  (n+l)C  p  V 
n  co 


2  2 

,*ru 
oo  is 


t  1/n+l 

XQ 


(2) 


where  the  integral  Iq  is  given  by 

s 

'q  *  ft?/?**1  d’  (3) 

-o 

Vaglio- Laurin,  R.,  11  Laminar  Heat  Transfer  on  Three-Dimensional  Blunt 
Nosed  Bodies  in  Hypersonic  Flow,  "  ARS  Journal,  29(2),  pp.  123-129 
(February  1959). 

7 

Vaglio- Laurin,  R.,  "Turbulent  Heat  Transfer  on  Blunt  Nosed  Bodies  in  Two- 
Dimensional  and  General  Three-Dimensional  Hypersonic  Flow,  "  Journal  of 
the  Aero/Space  Sciences,  27(1),  pp.  27-36  (January  I960). 

g 

Allen,  H,  Julian,  and  A.  J.  Eggers,  Jr.,  "A  Study  of  the  Motion  and  Aero¬ 
dynamic  Heating  of  Ballistic  Missiles  Entering  the  Earth's  Atmosphere  at 
High  Supersonic  Speeds,  "  NACA  Report  1381  (1958). 
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Allen  and  Eggers  integrated  an  expression  analogous  to  Eq.  (2),  representing 

Pand  V  as  functions  of  altitude  y  from  their  trajectory  analysis  and 
oo  oo 

assuming  ReQ  to  be  constant.  In  the  present  analysis,  the  shock  layer 

Reynolds  number  Re  is  represented  as  a  function  of  y.  Then,  following 

°  9 

the  general  integration  procedure  of  Brunner  ,  the  total  integrated  heat 
transfer  to  a  body  throughout  a  trajectory  is  given  for  laminar  flow  by 


QL=C^4«[^]‘/2erf(^) 


(4) 


and  for  turbulent  flow  by 


C 


T  T  9/5 
*E  B 


[BT]4'5 


/ 


-t 


5/4 


dt 


(5) 


The  complete  derivation  of  Eqs.  (4)  and  (5)  is  given  in  Ref.  10.  The  /? 
and  r _  in  these  equations  are  the  ballistic  coefficient  and  the  base  radius, 
respectively.  In  the  derivation  of  these  equations,  the  ballistic  coefficient 
has  been  assumed  to  be  constant  throughout  the  trajectory.  This  is  a  reason¬ 
ably  good  assumption  over  the  altitude  range  of  importance  for  convective 
heat  transfer,  i.e.,  less  than  100  kft.  Optimum  shapes  allowing  for  the 
changes  in  through  a  trajectory  have  been  reported  recently.  The 


9 

Brunner,  M.  J.,  "Analysis  of  Aerodynamic  Heating  for  a  Re-entrant  Space 
Vehicle,  "  Journal  of  Heat  Transfer,  8_1(8),  pp.  223-229  (August  1959). 

^Baker,  R.  L.,  and  R.F.  Kramer,  "Evaluation  of  Total  Heat  Transfer  in 
Hypersonic  Flow  Environments,"  Report  No.  TR-0077(2550- 15) -6,  The 
Aerospace  Corporation,  El  Segundo,  California. 

^Yelmgren,  Capt.  Kevin  E.,  "Optimum  Nose  Shape  for  Transpiration- 
Cooled  Reentry  Vehicles,  "  AFFDL  TR -78- 11 17,  Air  Force  Flight 
Dynamics  Laboratory  (November  1978). 
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other  parameters  are  given  by 


C  =  1.866x  10' 


(—) 

1/2 

L 

(  y+1) 

6 

\y- 1/ 

_T  4.70x10 

\Y-ll 

(2W+3)X 

'“T75- 

(2w+9) 

»lnvl12 


,L  .  E  E 

E  ~  .  .  *  ,1/2 
(sin0£) 


»l/5  9/5 

E  E 


T  4/5 


E  ,  .  -  .1/5 

(sin^E) 


(6a,  b) 


(7a,  b) 


'q  =  / p«u 


_ _ 2  ~  T 

p  r  ds,  ii 

ee  ’  Q 


5 

*  f  W 

•'a 


-  77  -5/4 
"  "  r  ds 
e 


(8a,  b) 


B  = 


P° 

XfisinOr 


B 


L  =  (2q>+3 )  D  bT 


(2&H-9) 

10 


B  (9a,  b,  c) 


The  variables  in  these  equations  are  the  air  specific  heat  ratio  and 
viscosity  law  exponent,  Y  and  (J;  the  density  scale  height  factor  in  the  expo¬ 
nential  atmosphere,  X;  the  reentry  velocity  and  entry  flight  path  angle,  VE 

and  0-,  the  air  viscosity  evaluated  at  the  entry  stagnation  temperature,  fj.-, 

^  L  T 

the  ballistic  factor  B  and  the  heat  transfer  integrals  Iq  and  Iq  which  depend 

upon  an  integration  in  body  surface  running  length  of  the  local  boundary  layer 

edge  density,  p g,  velocity  ug,  viscosity,  and  local  radial  coordinate  r. 

As  indicated  earlier,  previous  investigators  have  sought  to  minimize 
the  integrated  heat  transfer  to  a  body  for  the  case  in  which  the  freestream 
environment  is  constant.  Using  the  Vaglio-Laurin  convective  heat  transfer 
relationships  and  integrating  over  the  surface  area  as  in  Eq.  (1)  results  in 
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the  following  expression  for  the  integrated  heat  transfer  rate  to  a  body  in  a 

r  ,  .  ,1  10 

constant  environment  for  laminar  flow 


_ dQL/jt  =  1.328  (£l)  [Xq]  Reo1/2  (10a) 

p  V  H  (TTrl, ) 

rco  oo  o  \  B/ 

and  for  turbulent  flow 

,dQ.T/dt  =  0>072o(£t)  [Xq]  Reol/5  (10b) 

p  v  h  (jrrr.) 

oo  oo  o  \  B/ 


where  ReQ  is  a  shock  layer  Reynolds  number  defined  by 


Re  = 
o 


Pq(2H  )1/2rB 


(11) 


It  is  now  evident  that  if  we  fix  fg,  r&,  fi  and  B  in  Eq.  (4)  or  P^V ^ 

H  r  and  Re  in  Eq.  (10a),  then  the  total  heat  transfer  in  the  first  case 
o’  B  °  t  .  , 

and  the  total  heat  transfer  rate  dQ  /dt  in  the  second  case  will  be  minimized 

T  T 

if  we  minimize  the  integral  LT.  Similarly,  in  the  turbulent  case  Q  and 

'T  U  'p 

dQ  /dt  are  minimized  by  minimizing  the  integral  Iq.  The  relationship  of 
the  ballistic  trajectory  environment  problem  considered  in  this  work  to  the 
constant  environment  problem  considered  by  other  investigators  is  thus 
established. 

B.  EVALUATION  OF  iQand  Iq 

The  heat  transfer  integrals  Iq  and  Iq,  given  by  Eqs.  (8a)  and  (8b), 
respectively,  depend  intimately  upon  the  nose  tip  geometry,  i.  e.,  upon  the 
shape  contour  of  the  nose  tip.  As  shown  in  Fig.  1,  the  local  nose  tip  coordi¬ 
nates  r  and  z  and  the  surface  running  length  s  are  related  to  the  local  nose 
tip  surface  angle  a  by  the  following  simple  trigonometric  relationships 

=  tana,  4^  =  sinOI  (12) 

dz  ’  ds 
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To  a  reasonably  good  approximation  the  local  nondimens ional  boundary 

layer  edge  properties  p^,  and  in  Eqs.  (8a)  and  (8b)  may  be  represented 

as  functions  of  a.  In  the  present  work  the  simplest  equations  of  this  type 

have  been  used,  i.  e.,  the  hypersonic  Newtonian  approximations  for  both 

pressure  and  velocity.  Thus,  the  local  pressure  p  =  p  /p  is  given  by 

e  e  t2 


.  2 

sin  a 


(13) 


and  the  local  velocity  u^  =  n^/ yj  2Hq  is  given  by 


u  =  co  set 
e 


(14) 


If  the  viscosity  is  assumed  to  be  proportional  to  the  temperature  and  inde¬ 
pendent  of  pressure,  the  density- viscosity  product  is  then  directly  propor¬ 
tional  to  the  pressure,  i.  e. 


p  U  , 

—  —  e  e  . 

PH  =  ' —  =  sin  a 

e'e  p  p 
o  o 


(15) 


The  utility  rendered  by  the  simplicity  of  these  equations  was  felt  to 
compensate  for  inherent  inaccuracies  associated  with  them  for  purposes  of 
comparing  relative  heat  transfer  rates  to  various  nose  tip  shapes  as  discussed 
herein.  Work  is  currently  in  progress  to  reconsider  the  problem  discussed 
in  this  report  and  to  replace  the  approximations  given  by  Eqs.  (13)  through 
(15)  with  numerical  inviscid  flow  field  calculation  procedures.  Numerical 
flow  field  calculations  for  fixed  body  shapes  somewhat  similar  to  those 
discussed  here,  i.e.,  flat-face  cones,  are  given  in  Ref.  10. 


-12- 


c. 


LOCAL  PROPERTIES  ON  A  FLAT  FACE 


In  order  to  allow  for  solutions  which  may  have  a  flat  'ace,  it  is  neces¬ 
sary  to  have  an  alternate  expression  for  the  boundary  layer  edge  velocity  ug. 

This  is  because  direct  use  of  Eq.  (14)  would  indicate  that  u  and,  therefore, 

L  T  e 

IQ  and  Iq  are  all  identically  zero  on  a  flat  face.  In  Refs.  2  and  5,  the  follow¬ 
ing  approximation  for  the  product,  obtained  from  a  correlation  of 

method  of  integral  relations  solutions,  was  used. 


pun  — 

r  e  e^e  _  _r_ 

jU  —  _ **, 

t'VV'e1  r' 


(16) 


Equation  (16)  indicates  that  for  a  flat-faced  body,  the  normalized  Peue^e 
product  varies  linearly  with  distance  to  the  sonic  point.  The  sonic  point 
values  of  the  dependent  variables  may  be  obtained  by  isentropic  expansion 
from  the  stagnation  conditions.  The  terms  p  ,  u  and  p  in  Eq.  (16)  may 

then  be  renormalized  by  the  stagnation  values  Pq,  J 2Hq,  Using  ideal 

gas  relationships,  we  obtain 


p  u  u 
re  e  e 

',o(2Ho)1/2"o 


(17a) 


F  (Y) 


(17b) 


In  Section  III,  Eqs.  (12)  through  (17)  are  used  to  express  the  heat  transfer 
integral  Iq  explicitly  in  terms  of  the  body  geometry. 
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1 


Since  T  is  defined  to  be  the  ratio  of  the  diameter  to  the  length  of  the  body, 

y(2/r)  =  1  (22) 

We  wish  to  determine  the  function  y(x),  chosen  from  a  suitable  class  of 
functions  on  the  interval  0  <  x  <  2/r  which  minimizes  (or  maximizes)  the 
functional  Iq,  Eq.  (21),  and  which  satisfies  the  condition  (22).  Any  such 
minimum  (or  maximum)  must  satisfy  the  Euler  equation  associated  with  this 
problem. 


yy(l-3y2] 


(n+1)  .2  , .  .4, 

t1_y  J 


(23) 


The  boundary  condition  y(2 /t)  =  1  is  given,  but  y(0)  is  not  specified.  Hence, 
this  is  a  variable  end  point  problem  of  the  calculus  of  variations,  and  the 
following  transversality  condition  must  be  satisfied  at  x  =  0. 


n+1 


<i+V 


n+2 
2  (n+3) 


F(y) 


■  0 


(24a) 


Since  the  integrand  of  Eq.  (21)  is  not  explicitly  a  function  of  x,  a  first 
integral  of  the  Euler  equation  (23)  can  be  derived.  This  equation  may  be 
written  in  the  following  form: 


n+1  C(l+y2) 

y  7  7 

y  (y  -i) 


(24b) 


Examination  of  this  equation  reveals  that  there  is  no  real  value  of  y  which 
allows  y  to  be  zero.  Hence,  there  is  no  solution  of  the  Euler  equation  for 
which  y(0)  *  0. 
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If  a  flat-faced,  y(0)  *  0,  solution  exists,  then  [from  Eq.  (24a)]  the  initial 
slope  must  be  a  solution  of  the  equation 


y(Q) 

[i  +y(0)2]2 


_  (n+2) 

~  2(n+3) 


F(y) 


(25) 


This  equation  is  a  quartic  in  y(0).  An  examination  of  this  quartic  for  the 
range  of  parameters  (n,y)  of  interest  reveals  two  real  positive  roots  and  a 
pair  of  complex  conjugate  roots.  Thus,  there  are  apparently  two  functions 
y(x)  which  satisfy  the  Euler  equation  and  the  trans versality  condition.  It 
should  be  observed  that  the  values  of  the  two  permissible  initial  slopes  are 
independent  of  r. 

One  might  expect  that  the  two  solutions  for  y(x)  discussed  in  the 
preceding  paragraph  represent  a  minimum  and  a  maximum  of  Iq.  There 
are  further  checks  which  may  be  made  to  verify  or  limit  this  conclusion. 

The  necessary  condition  of  Legendre  implies  that  one  of  the  two  solutions 
for  y(x)  is  a  minimum  (maximum)  if  y(x)  <  n/3/3  [y(x)  >  ”3/3]  for  all  x  in  the 
interval  0  <  x  <  2 /t.  The  Jacobi  condition  is  somewhat  difficult  to  apply. 
However,  it  can  be  shown  that  the  necessary  condition  of  Weierstrass  cannot 
be  satisfied  over  the  class  of  all  possible  shapes.  Hence,  neither  of  the 
above  solutions  corresponds  to  a  minimum  or  a  maximum  of  Iq  over  the 
class  of  all  piecewise  continuously  differentiable  functions  on  the  interval 
0  <  x  <  2/T  which  satisfy  the  boundary  condition  Eq.  (22).  Hull^  has  pointed 
out  a  similar  occurrence  in  a  related  problem. 

B.  INTERPRETATION  OF  SOLUTIONS 

The  solutions  discussed  in  the  preceding  subsection  are  insufficient  as 
global  extrema.  This  may  be  demonstrated  as  follows:  The  integral  appearing 
in  Eq.  (19)  is  zero  over  any  segment  for  which  y  is  zero  or  for  which  y  is 
infinite.  Hence  as  y(0)  approaches  zero,  there  are  infinite  numbers  of  shapes 
for  which  Iq  is  arbitrarily  close  to  zero,  i.e.,  shapes  composed  of  vertical 
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and  horizontal  segments.  Further  consider  functions  y(x)  which  are  straight 
lines  between  y(0)  at  x  =  0  and  1  at  x  =  2 /t.  As  y(0)  becomes  arbitrarily 
large,  Iq  will  become  arbitrarily  large. 

The  two  solutions  which  satisfy  the  Euler  equation  [Eq.  (23)],  the 
transversality  condition  [Eq.  (25)]  and  the  boundary  condition  [Eq.  (22)] 
will  be  referred  to  as  a  "minimum"  solution  or  a  "maximum"  solution  pro¬ 
vided  they  satisfy  the  Legendre  condition  for  a  maximum  or  a  minimum.  As 
has  been  shown,  these  solutions  do  not  represent  minima  or  maxima  over 
all  possible  shapes.  However,  a  minimum  solution  (y$  *s/3/3)  will  satisfy 
the  Weierstrass  condition  for  a  minimum  over  the  class  of  all  piecewise 
continuously  differentiable  functions  y(x)  on  the  interval  0<x$2/r  for 
which  y(2/T)  =  1  and  for  which  y(x)  <  \/3/3  for  all  x  in  the  interval.  Similarly, 

a  maximum  solution  (y  >  ^/3)  will  satisfy  the  Weierstrass  condition  for  a 
maximum  over  the  class  of  all  piecewise  continously  differentiable  functions 
y(x)  on  the  interval  0  <  xS  2/r  for  which  y(2 /t)  =  1  and  for  which  y(x)  >  */3/3 
for  all  x  in  the  interval.  With  these  restrictions  on  derivatives,  not  only  is 
the  Weierstrass  condition  satisfied,  but  the  Weierstrass -Erdmann  corner 
condition  cannot  be  satisfied  by  any  pair  of  admissible  slopes.  Hence,  there 
are  no  minimum  or  maximum  solutions  with  discontinuous  derivatives. 

These  limits  seemed  unduly  restrictive.  However,  after  obtaining  the 
solutions  it  was  found  that  they  could  be  extended.  For  instance,  it  was 
observed  that  the  calculated  maximum  solutions  had  the  property  that  y(x)  £  1 
for  all  x.  Examination  of  the  Weierstrass  excess  function  for  this  problem 
leads  to  the  following  conclusion:  If  a  maximum  solution  is  found  whose 
derivative  is  greater  than  or  equal  to  one  for  all  x,  then  the  Iq  corresponding 
to  this  solution  will  be  larger  than  the  Iq  corresponding  to  any  piecewise 
continuously  differentiable  function  y(x)  on  the  interval  0  <  x  5  2/r  for  which 
y (2/r)  =  1  and  for  which  y(x)  >  0  for  all  x  in  the  interval. 
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When  a  minimum  solution  is  found,  its  maximum  derivative  y  over 

niftx 

the  interval  0  £  x  <  2/T  may  be  determined.  If  the  constant  A  is  defined  by 


A 


1 


•  2 
y 

7  max 


2y. 


max 


(26) 


then  examination  of  the  Weierstrass  E-function  for  the  problem  leads  to  the 
following  situation.  If  a  minimum  solution  is  found  with  maximum  derivative 
v  and  constant  A,  then  the  1,^  corresponding  to  this  solution  will  be 
smaller  than  the  Iq  corresponding  to  any  piecewise  continuously  differentiable 
function  y(x)  on  the  interval  0  ^  x  £  2/T  for  which  y(2 /T)  =  1  and  for  which 
y(x)  <  A  for  all  x  in  the  interval.  It  is  in  this  limited  sense  that  the  solutions 
given  in  the  following  sections  may  be  considered  minima  and  maxima. 


C  COMPUTATION  OF  SOLUTIONS 

The  actual  computation  of  a  solution  consists  of  solving  a  two-point 
boundary  value  problem  for  a  second-order  nonlinear  ordinary  differential 
equation.  The  differential  equation  is  the  Euler  equation  (23).  The 
boundary  conditions  are  y(2/T)  =  1  and  at  x  =  0  the  trans versality  condition 
(25).  A  program  was  written  to  solve  a  class  of  such  problems,  including 
the  one  discussed  herein.  In  this  program  the  differential  equations  and 
boundary  conditions  were  replaced  by  finite  difference  approximations,  and 
the  resulting  equations  were  solved  by  Newton's  method  (quasilinearization). 
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IV.  RESULTS  AND  DISCUSSION 


Results  obtained  by  numerical  solution  of  the  Euler  equation  with  ap¬ 
propriate  boundary  conditions,  as  discussed  above,  are  presented  here  for 
the  cases  of  both  laminar  and  turbulent  flow  in  the  boundary  layer,  i.  e., 
n  =  1  and  n  =  1/4,  respectively.  In  all  cases  the  specific  heat  ratio  Y  has 
been  assumed  to  be  the  hypersonic  value  of  1.2.  In  the  first  subsection  below, 
the  body  shapes  and  minimum  heat  transfer  integrals  Iq  and  ^  are  shown 
for  a  wide  range  of  the  fineness  ratio  T.  Then,  in  the  next  subsection,  the 
minimum  heat  transfer  results  are  compared  with  heat  transfer  rates  to 
other  selected  body  shapes  and  to  some  maximum  heat  transfer  solutions. 

A.  MINIMUM  HEAT  TRANSFER  SOLUTIONS 

Laminar  flow  minimum  heat  transfer  shapes  for  several  values  of  the 
fineness  ratio  T  are  shown  in  Fig.  2  where  the  shapes  as  they  appear  physi¬ 
cally  are  shown  in  the  upper  part  of  the  figure.  In  the  lower  part  of  the 
figure,  the  shapes  are  shown  with  the  radial  coordinate  normalized  by  the 
base  radius  and  the  axial  coordinate  normalized  by  the  body  length.  In  this 
way  the  maximum  value  of  both  r  =  r/r^  and  z  =  z/£  is  unity  since  two 
reference  lengths  have  been  used.  In  the  remainder  of  this  report,  all  body 
shapes  will  be  presented  normalized  as  in  the  lower  part  of  Fig.  2.  It  should 
be  kept  in  mind  that  presenting  results  in  this  way  masks  the  needlelike 
appearance  of  shapes  with  very  small  fineness  ratios  T  as  well  as  the  blunt 
disklike  appearance  of  shapes  with  large  fineness  ratios. 

All  of  the  shapes  in  Fig.  2,  for  nonzero  values  of  T,  exhibit  a  flat  face. 
The  ratio  of  the  radius  of  the  flat  face  to  that  of  the  base  increases  with 
increasing  T.  In  the  limit  of  T  approaching  infinity,  the  shape  becomes  a 
flat-faced  cylinder  of  zero  length,  i.e.,  a  disk.  In  the  opposite  limit,  as  T 
approaches  zero,  the  flat-face  height  approaches  zero,  and  the  shape 
approaches  the  shape  obtained  invoking  the  slender  body  approximation,  i.e., 
a  1/2  power  law  body  .  However,  since  the  slender  body  power  law  shape  is 


only  approached  for  fineness  ratios  less  than  0.02,  this  shape  appears 
to  be  of  little  practical  interest. 

In  Fig.  3,  turbulent  flow  minimum  heat  transfer  shapes  are  shown  and 

compared  with  the  laminar  shapes  for  corresponding  fineness  ratios.  As 

5 

concluded  by  others  ,  the  turbulent  shapes  are  surprisingly  similar  to 

the  laminar  shapes  for  a  given  fineness  ratio.  In  the  limit  of  large  T,  the 

turbulent  minimum  heat  transfer  shape  also  approaches  a  flat-faced 

cylinder;  in  the  limit  of  T  approaching  zero  it  can  be  easily  shown  that  the 

turbulent  minimum  heat  transfer  shape  is  an  8/13  power  law  body.  The 

slender  body  derivation  of  power  law  shapes  for  minimum  laminar  and 

turbulent  heat  transfer  is  given  in  Appendix  A.  Thus,  in  the  slender  body 

limit  the  minimum  laminar  heat  transfer  shape  is  a  1/2  power  law  body, 

the  minimum  turbulent  heat  transfer  shape  is  an  8/13  power  law  body,  and 

(from  Ref.  12)  the  minimum  drag  shape  is  a  3/4  power  law  body.  The  non- 

slender  body  minimum  heat  transfer  shapes  found  in  the  present  work 

approach  the  slender  body  minimum  heat  transfer  shapes  only  at  very  small 

fineness  ratios  (T  <  0.02).  However,  the  nonslender  body  minimum  drag 

1 2 

shapes  approach  a  3/4  power  law  body  at  much  larger  fineness  ratios 
(TWO. 2). 

The  question  of  how  much  reduction  in  the  heat  transfer  is  achieved  by 
minimum  heat  transfer  shape  is  considered  in  the  results  presented  in 
Figs.  4  through  6  and  in  the  following  subsections.  Recall  that  Eqs,  (4) 
and  (5)  represent  expressions  for  the  integrated  heat  transfer  to  a  nose  tip 
throughout  a  trajectory.  Likewise,  Eqs.  (10a)  and  (10b)  represent  similar 
expressions  for  the  constant  freestream  environment  case.  For  all  cases, 
trajectory  environment  or  constant  freestream  environment,  the  parameter 
needed  to  obtain  numerical  predictions  of  heat  transfer  is  the  heat  transfer 
integral  or  body  shape  parameter  Iq. 

__ 

Eggers,  A.J.,  Jr.,  M.  M.  Resnikoff  and  D.  H.  Dennis,  "Bodies  of  Revolu¬ 
tion  Having  Minimum  Drag  at  High  Supersonic  Speeds,  "  NACA  Report 
No.  1306  (1957). 
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Minimum  I~  as  a  Function  of  Fineness  Ratio 


b)  TURBULENT 


In  Fig.  4,  IQ  for  the  minimum  laminar  and  turbulent  heat  transfer 

shapes  is  shown  as  a  function  of  T.  The  Iq  for  both  laminar  and  turbulent 

flow  is  quite  insensitive  to  fineness  ratio  for  1  <  T  <  10.  This  is  because 

the  ratio  of  flat-face  height  to  base  radius  is  greater  than  0.  9  for  all  T>  1. 

Thus,  for  T  >1,  a  large  fraction  of  the  total  heat  transfer  occurs  over  the 

L  T 

flat  face.  However,  for  T  <  1  the  Iq  and  Iq  values  decrease  significantly  as 
the  body  or  nose  tip  becomes  more  slender,  i.e.,  decreasing  T.  This 
decrease  is  due  in  large  part  to  a  significant  reduction  in  the  flat-face  height, 
which  decreases  the  area  associated  with  the  relatively  high  flat-face  heat 
transfer. 

L  T 

In  the  limit  of  T approaching  zero,  I0  and  I~  both  go  to  zero.  However, 

Li  T  w  w 

Iq/T  and  Iq/t  are  finite  for  T  =  0,  and  the  limiting  values  are  given  by  slen¬ 
der  body  theory  (see  Appendix  B).  They  are 


I  Ix 
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(27a,  b) 


t  t 

Comparing  the  numerical  values  of  Iq  and  Iq  from  these  equations  with  the 
nons lender  body  calculated  values  in  Fig.  4,  we  see  that  for  T$.  0.  02  the 
slender  and  nonslender  body  calculated  values  are  virtually  identical.  A 
quantitative  comparison  of  minimum  heat  transfer  rates  to  those  for  other 
body  shapes  is  given  below. 

B.  COMPARISON  WITH  OTHER  BODY  SHAPES 

The  minimum  heat  transfer  body  shapes  and  heat  transfer  integrals 
discussed  above  will  now  be  compared  with  these  quantities  for  a  sphere, 
truncated  spheres,  sharp  cones,  minimum  drag  bodies,  stable  ablating 
shapes  and  maximum  heat  transfer  bodies.  Except  for  the  spheres,  sharp 
cones  and  minimum  drag  bodies,  the  remaining  body  shapes  listed  above 
depend  upon  the  boundary  layer  state,  i.e.,  laminar  or  turbulent. 
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The  minimum  and  maximum  heat  transfer  solutions  were  obtained 
directly  by  numerical  solution  of  the  Euler  equation  as  discussed  previously. 
The  minimum  drag  body  solutions  were  obtained  by  application  of  the  present 
numerical  calculation  procedure  to  the  appropriate  equations  as  given  in 
Ref.  12.  A  stable  ablating  shape  is  one  for  which,  even  though  it  is  ablating 
and  its  surface  is  receding,  the  body  shape  contour  is  independent  of 
time  The  following  analytic  equations  for  laminar  and  turbulent  stable 
ablating  nose  tip  shapes,  using  the  approximations  given  by  Eqs.  (13)  through 
(15),  are  derived  in  Appendix  C.  For  laminar  flow 


/  1-r2  +  1 


(28a) 


and  for  turbulent  flow 
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The  shape  contours  for  the  six  body  types  listed  above  are  compared  in 
Figs.  5a  and  5b  for  fineness  ratios  of  8  and  4  for  laminar  and  turbulent  flow, 
respectively.  The  laminar  stable  ablating  nose  tip  shape  in  Fig.  5a  most 
closely  resembles  the  truncated  sphere  shape,  whereas  the  turbulent  stable 
shape  in  Fig.  5b  is  similar  to  the  maximum  heat  transfer  shape.  The  flat- 
nose  height  is  smallest  for  the  maximum  heat  transfer  body,  somewhat 
larger  for  the  minimum  drag  body  and  largest  for  the  minimum  heat  transfer 
body. 


Baker,  R.  L.,  "The  Effect  of  Freestream  and  Material  Parameters  on 
Stable  Ablating  Nosetip  Shapes,  n  Report  No.  TOR-0074(4450-76)-10,  The 
Aerospace  Corporation,  El  Segundo,  California  (May  1974). 
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The  heat  transfer  integrals  1q  and  Iq  are  shown  for  these  body  shapes 
in  Figs.  6a  and  6b,  respectively.  The  relative  Iq  values  for  the  different 
body  shapes  are  generally  similar  for  laminar  and  turbulent  flow.  The  maxi¬ 
mum  solution  Iq  values  are  up  to  four  times  greater  than  the  minimum  solu¬ 
tion  values  with  this  ratio  decreasing  as  T  increases  for  2  <  T<  8.  This 
means  that  for  this  range  of  fineness  ratio,  the  heat  transfer  for  the  laminar 
maximum  solutions  is  up  to  two  times  greater  than  the  laminar  minimum 
solution  values.  The  corresponding  heat  transfer  increase  for  the  minimum 
to  maximum  turbulent  solutions  is  up  to  three  times  greater. 

The  Iq  values  for  truncated  spheres,  sharp  cones,  stable  ablating  shapes 
and  minimum  drag  bodies  all  lie  between  the  minimum  and  maximum  solution 
values  but  tend  to  be  much  closer  to  the  maximum  value  for  a  given  T.  An 
exception  is  the  Iq  for  minimum  drag  bodies  for  small  T.  The  Iq  for  these 
bodies  peaks  for  T  =  2  and  decreases  for  all  other  fineness  ratios.  At  T  *  2, 
the  minimum  drag  body  Iq's  for  laminar  and  turbulent  flow  are  close  to  the 
corresponding  maximum  heat  transfer  solution  values.  However,  for  T 
decreasing  and  approaching  zero  (T  <  0.  02),  the  minimum  drag  body  Iq  values 
approach  the  minimum  heat  transfer  body  values.  Maximum  heat  transfer 
solutions  could  not  be  found  for  T  <  2. 

C.  SOLUTION  CHARACTERISTICS  AND  ACCURACY 

Since  many  of  the  extremal  solution  body  shapes  discussed  here  have  a 
flat  nose  followed  by  a  body  contour,  it  is  of  interest  to  compare  these  solu¬ 
tions  in  terms  of  the  body  angle  immediately  aft  of  the  flat  face  and  the  ratio 
of  flat- face  height  to  base  radius.  For  minimum  heat  transfer,  minimum 
drag  and  maximum  heat  transfer  bodies,  the  body  angle  aft  of  the  flat  face  is 
about  4,  45,  and  63  deg,  respectively,  independent  of  T  and  essentially  inde¬ 
pendent  of  the  boundary  layer  state,  i.e.,  laminar  or  turbulent.  This  is  shown 
in  Table  1,  Caution  must  be  exercised  concerning  the  maximum  solutions. 
While  the  local  flow  at  the  expansion  aft  of  the  flat  face  is  indicated  to  be 
supersonic  using  the  present  approximations,  more  exact  calculations  could 


indicate  subsonic  flow  at  this  location.  The  present  formulation  assumes  the 
flow  to  be  supersonic  at  the  expansion  corner.  If  this  condition  is  not  met, 
the  proof  of  the  existence  of  maximum  solutions  becomes  much  more  difficult 
and  the  present  solutions  are  not  valid. 

Table  1.  Body  Angle  Aft  of  Expansion  Corner  for 
Different  Variational  Calculus  Solutions 


Laminar 

Turbulent 

Minimum  Heat  Transfer 
Solution 

3°  41' 

3° 

24' 

Minimum  Drag  Solution 

45° 

45° 

Maximum  Heat  Transfer 

65°  39' 

66° 

23' 

The  flat  face-base  radius  ratio  for  minimum  heat  transfer,  maximum 
heat  transfer  and  minimum  drag  shapes  is  shown  in  Fig.  7  a  function  of 
fineness  ratio.  This  ratio  increases  from  essentially  zero  to  0.6-0.  7,  for  T 
increasing  from  2  to  10  for  the  maximum  heat  transfer  solutions.  For  the 
minimum  drag  body,  the  flat-face  height  increases  over  the  same  range  for 
fineness  ratios  from  0.  2  to  10,  and  for  the  minimum  heat  transfer  body 
the  T  range  is  about  0.02  to  10  for  the  same  increase  in  flat-face  height. 

A  final  word  of  caution  is  in  order  regarding  the  accuracy  of  the  pre¬ 
dicted  reduction  in  heat  transfer  associated  with  the  minimum  heat  transfer 
shapes.  Calculations  using  numerical  inviscid  flowfield  computer  codes  for 
flat-face  cones  were  carried  out  recently  by  Baker  and  Kramer1  .  These 
calculations  indicate  that  while  the  local  flow  along  the  surface  is  highly  over- 
expanded  immediately  aft  of  the  flat-face  expansion  corner,  the  flow  rapidly 
recompresses  to  pressure  levels  considerably  above  those  for  the  same  sta¬ 
tion  on  a  sphere-cone.  Because  of  this  behavior,  the  reduction  in  heat 
transfer  attributable  to  minimum  heat  transfer  type  shapes,  when  calculated 
using  numerical  flowfield  calculation  methods,  is  expected  to  be  a  factor  of 
up  to  2  rather  than  the  factor  of  3  predicted  by  the  present  work. 
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V.  SUMMARY  AND  CONCLUSIONS 


Analytic  expressions  for  total  integrated  heat  transfer  to  a  nose  tip 
through  a  trajectory  and  variational  calculus  procedures  have  been  used  to 
determine  nose  tip  shape  contours  having  minimum  total  heat  transfer  for 
specified  fineness  ratio.  Solutions  for  both  laminar  and  turbulent  boundary 
layer  flow  were  obtained  by  numerical  solution  of  the  appropriate  Euler  equa¬ 
tion  and  boundary  conditions. 

For  fineness  ratio  T  greater  than  zero,  all  solutions  found  have  a  flat 
face.  The  ratio  of  flat-face  height  to  the  base  radius  increases  as  T  increases. 
As  T  becomes  large  (T  >10),  both  the  laminar  and  turbulent  solutions  approach 
a  flat-faced  cylinder  of  vanishing  length,  i.e.,  a  disk.  As  T  approaches  zero 
(T  <  0.02),  the  solutions  approach  the  slender  body  (T  =  0)  solutions  which 
are  a  1/2  power  law  body  for  laminar  flow  and  an  8/13  power  law  body  for 
turbulent  flow.  For  a  constant  base  radius,  the  relative  heat  transfer  rates 
to  a  family  of  minimum  heat  transfer  shapes  decreases  monotonically  as  T 
decreases,  i.e.,  as  the  bodies  become  more  slender.  For  T  > 0. 02,  the  lami¬ 
nar  and  turbulent  shapes  are  surprisingly  similar  to  one  another,  i.e.,  the 
flat- face  height  to  base  radius  ratio  and  the  afterbody  shape  contour  are 
essentially  independent  of  the  boundary  layer  state. 

Additional  solutions  to  the  Euler  equation  depicting  maximum  heat 
transfer  body  shapes  were  found  for  T  >  2.  The  heat  transfer  to  the  maximum 
solution  shapes  is  up  to  two  times  greater  than  for  the  minimum  solution 
shapes  for  laminar  flow  and  up  to  three  times  greater  for  turbulent  flow  for 
2  <  T<  8.  Laminar  and  turbulent  heat  transfer  rates  to  sharp  cones,  spheres, 
truncated  spheres  and  stable  ablating  shapes  were  found  to  be  generally  much 
closer  to  the  maximum  solution  values  than  to  the  minimum  solution  values. 
The  laminar  and  turbulent  heat  transfer  to  a  minimum  drag  body  was  also 
calculated. 
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The  minimum  and  maximum  solutions  obtained  could  be  shown  mathe¬ 
matically  to  be  minima  and  maxima  only  over  restricted  classes  of  body 
shapes.  However,  of  all  the  physically  reasonable  additional  body  shapes 
considered,  none  were  found  which  result  in  heat  transfer  rates  lower  than 
the  minimum  solution  or  greater  than  the  maximum  solution  values. 
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APPENDIX  A 


SLENDER  BODY  MINIMUM  HEAT  TRANSFER  SHAPES 


Following  the  methods  and  nomenclature  of  Ref.  1,  define  the  following 
nondimensional  coordinates: 


T]  =  Y,  t  =  7  x 


2R, 


r  = 


(A  - 1) 


Then,  after  differentiation  and  direct  substitution  into  Eq.  (17),  we  have 
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Since  T  «  1  for  a  slender  body,  in  the  slender  body  approximation  the 
expression  for  Iq  becomes 
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(Slender  Body)  IQ  =  ~  J  r}n  +  1  T]Zd$  =  J  F{rj,  £)d£  (A-3) 


For  the  case  in  which  the  function  F  is  not  an  explicit  function  of  the 
independent  variable,  | ,  the  first  integral  of  the  Euler  equation  can  be  written 
immediately  as 
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Performing  the  indicated  differential  and  algebraic  manipulations,  the  follow¬ 
ing  first  integral  of  the  Euler  equation  for  this  problem  is  obtained  from 
Eqs.  (A-3)  and  (A-4)  after  rearrangement 


c„-,n+I)/2 


The  boundary  conditions  for  Eq.  (A -5)  are 


(A -5) 


>7(0)  =  0 

>7(D  =  1  (A -6) 


Integrating  Eq.  (A-5)  and  applying  the  boundary  conditions,  we  obtain 
for  laminar  flow  (n=l) 


T]  * 


(A -7) 


and  for  turbulent  flow  (n=l/4) 
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APPENDIX  B 

HEAT  TRANSFER  INTEGRA LS- -SLENDER  BODY  THEORY 

L  T 

The  values  of  the  heat  transfer  integrals  Iq  and  Iq  both  go  to  zero  as 
the  fineness  ratio,  T,  goes  to  zero.  However,  from  Eq.  (A -3),  we  have 
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For  laminar  flow  this  becomes 


Similarly,  for  turbulent  flow  Eq.  (B-l)  becomes 


(B-5) 


From  Eq.  (A -8),  for  turbulent  flow 
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Substituting  and  integrating,  we  have 


(B-7) 


Comparing  the  numerical  values  of  Iq  and  Iq  from  Eqs.  (B-4)  and 
(B-7)  with  the  non-slender  body  calculated  values  in  Fig.  5,  we  see  that  for 
T<  0.  02  the  slender  and  non-slender  body  calculated  values  are  very  close  to 
one  another. 


APPENDIX  C 


STABLE  ABLATING  NOSETIP  SHAPES 


The  transient  shape  change  equation  derived  in  Ref.  10  becomes  after 
differentiation  with  respect  to  s  and  appropriate  substitutions  related  to  the 
geometry  of  Fig.  1 
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For  a  steady-state  condition,  i.e.,  a  stable  ablating  nosetip  shape,  the  local 
body  angle.  Of,  is  not  a  function  of  time.  Thus 
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Integ  rating 


Q  sinCf 


=  constant  =  C 
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We  wish  to  determine  the  body  shape  such  that  the  heat  transfer  distri¬ 
bution  q  =  q(Ct,  s)  will  satisfy  Eq.  (C-3).  From  Ref.  10,  for  laminar  flow  in 
the  boundary  layer 
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Combining  Eqs.  (C-3)  and  (C-4)  and  assuming  that  the  effective  heat  of 

& 

ablation,  Q  ,  is  constant 
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Let 


f  = 


p  u  u 
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Then  Eq.  (C-5)  becomes 
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By  direct  substitution,  we  find  that  Eq.  (C-7)  is  satisfied  by 
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Combining  Eqs.  (C-6)  and  (C-8)  together  with  Eqs.  (11)  through  (13)  and  the 
trigonometric  relations 


QZ 


the  following  differential  equation  for  the  laminar  stable  ablating  nosetip 
shape  is  obtained; 
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where  A1j=  (C1QV/CL)  ,  or  defining  r  =  r/AJ 

1  -  >/l  -  r2 


-2 

r  = 


(C-lOb) 


Integrating  Eq.  (C-10)  and  applying  the  boundary  condition  z(0)  =  0,  we  obtain 
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For  turbulent  flow  in  the  boundary  layer  we  have  from  Ref.  11 
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Applying  the  same  procedures  as  in  the  laminar  case  above,  we  find  that  the 
differential  equation  describing  the  turbulent  stable  ablating  shape  is 
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Integrating  Eq,  (C-13)  with  the  aid  of  the  substitution  X  =  r  ,  applying  the 
boundary  condition  z(0)  =  0  and  nondimensionalizing,  we  obtain 
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